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ABSTRACT 

Previous  studies  have  shown  that  forcing  a  low  Reynolds 
number,  confined  2-D  wake  may  lead  to  greatly  enhanced 
molecular  mixing.  This  study  presents  whole-field 
measurements  of  the  velocity  and  vorticity  fields  in  such  a  flow 
in  order  to  examine  its  underlying  flow  structure  and  the 
features  connected  to  mixing  enhancement.  The  variation  of 
the  spanwise  vorticity  field  versus  forcing  amplitude  is 
presented.  Results  show  that  the  downstream  location  where 
the  magnitude  of  the  mean  velocity  gradient  5w/5z  starts  to 
increase  coincides  closely  with  the  location  where  mixing 
enhancement  has  been  noted  to  increase.  The  downstream 
evolution  of  the  streamwise  vorticity  shows  that  the  high 
amplitude  forcing  cases  corresponding  to  significant  mixing 
enhancement  are  characterized  by  large  values  of  co^.  For  these 
cases,  the  regions  of  move  quickly  away  from  the  facility 
sidewalls  towards  the  center-span  of  the  test  section.  In 
contrast,  these  regions  remain  close  to  the  sidewalls  in  the  low 
amplitude  forcing  conditions  which  result  in  a  modest  mixing 
increase. 


INTRODUCTION 

Forcing  a  low  Reynolds  number,  confined  2-D  wake  can  lead 
to  a  highly  three-dimensional  flow  and  a  large  increase  in 
mixing.  Estimates  of  the  amount  of  molecularly  mixed  fluid  in 
a  liquid-phase  flow  made  by  Koochesfahani  and  Nelson  (1997) 
have  shown  that  the  mixed  fluid  fraction  in  this  forced,  low 
Reynolds  number  flow  is  2.5  to  3  times  larger  than  that  seen  in 
high  Reynolds  number,  natural  two-stream  liquid  mixing  layers 
and  60%  larger  than  that  in  gas-phase  turbulent  shear  layers. 

Figure  1  shows  the  downstream  development  at  midspan  of 
the  mixed  fluid  fraction,  6p/6i,  at  a  fixed  frequency  for  three 
different  forcing  amplitudes  from  Nelson  and  Koochesfahani 
(1997).  The  mixed  fluid  fraction  is  defined  here  as  the  ratio  of 
product  thickness  measured  by  a  chemically  reacting  LIF 


technique  to  the  local  wake  width.  Note  in  Figure  1  that 
forcing  at  low  and  moderate  amplitudes  results  in  a  modest 
increase  in  the  amount  of  mixed  fluid.  The  high  amplirnde 
forcing  case,  however,  results  in  a  significant  increase  in  6p/6i. 
The  leveling  out  of  the  profile  for  x  >  20  cm  is  caused  by  the 
confinement  imposed  by  the  top  and  bottom  walls  of  the  flow 
facility. 

Earlier  studies  by  Roberts  (1985)  show  that  mixing 
enhancement  is  also  seen  when  the  forcing  frequency  is  varied. 
Significant  increases  in  the  amount  of  mixed  fluid  were  found 
in  this  study  when  the  flow  was  forced  at  twice  the  natural 
shedding  frequency.  The  mixing  studies  of  Roberts, 
Koochesfahani  and  Nelson,  and  MacKinnon  and 
Koochesfahani  (1997)  have  all  noted  that,  for  cases  of  large 
mixing  enhancement,  the  mixing  regions  start  near  the  side 
walls  of  the  facility  and  move  toward  center-span  farther 
downstream.  A  simple  vortex  model  was  proposed  by  Roberts 
to  explain  this  feature.  The  mixing  regions  near  the  side  walls 
were  modeled  by  a  pair  of  counter-rotating  streamwise  vortices 
formed  as  a  result  of  the  reorientation  of  the  primary  spanwise 
wake  vortices  due  to  the  side-wall  boundary  layer.  The 
induced  flow  field  of  this  vortex  system,  and  its  images,  cause 
the  lateral  movement  of  these  vortical  regions  toward  the 
center- span. 

All  previous  studies  of  the  wake  flow  discussed  here  have 
concentrated  on  either  flow  visualization  or  quantification  of 
the  mixing  characteristics.  While  all  these  studies  connect  the 
increase  in  mixing  to  the  downstream  evolution  of  the 
streamwise  vorticity,  no  information  has  been  available  on  the 
behavior  of  the  vorticity  field  in  this  flow.  In  the  present 
study,  the  vorticity  evolution  is  quantified  using  Molecular 
Tagging  Velocimetry  (MTV),  and  the  structure  of  the  velocity 
and  vorticity  fields  are  explored  in  the  context  of  previously 
observed  mixing  enhancement. 
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EXPERIMENTAL  FACILITY 

The  experiments  were  conducted  in  a  gravity  driven,  liquid 
two-stream  mixing  layer  facility  located  in  the  Turbulent 
Mixing  and  Unsteady  Aerodynamics  Laboratory  at  Michigan 
State  University.  This  is  the  same  facility  used  in  the  mixing 
studies  conducted  by  MacICinnon  and  Koochesfahani  (1997), 
Koochesfahani  and  Nelson  (1997),  and  Nelson  (1996).  The 
“high”,  “middle”,  and  “low”  amplitude  forcing  in  these  studies 
(see  also  Figure  1)  correspond  to  velocity  perturbation 
amplitudes  of  1 1.6%,  7.3%,  and  2.1%.  In  this  facility,  fluid 
pre-mixed  with  chemicals  necessary  for  the  molecular  tagging 
diagnostic  technique  are  pumped  into  two  constant  head 
reservoir  tanks  approximately  3.5  m  above  the  test  section. 
The  operational  parameters  of  the  experiments  were  selected 
to  be  the  same  as  previous  mixing  studies  in  order  to  allow 
direct  comparison.  The  free-stream  speeds  of  the  two  sides 
were  set  to  u^  =  Uj  =  9.4  cm/s  to  generate  a  wake  flow  with  an 
initial  Reynolds  number  of  about  100  based  on  the  wake 
momentum  thickness. 

Two-dimensional  perturbations  were  introduced  into  the 
upper  free-stream  of  this  flow  by  means  of  an  oscillating 
bellows  mechanism  in  the  upper  supply  line.  The  mechanism 
was  driven  by  an  electromagnetic  shaker  whose  command 
signal  originated  from  a  standard  function  generator.  The 
command  signal  was  recorded  for  each  experiment  to  facilitate 
the  phase-averaging  of  the  measured  data  which  will  be 
discussed  in  the  next  section.  Sinusoidal  perturbation  of  the 
flow  at  different  amplitudes  and  a  fixed  frequency  of  6 
which  is  approximately  the  natural  shedding  frequency,  were 
examined  in  this  study.  It  should  be  noted  that  even  though 
forcing  is  only  applied  to  the  upper  stream,  oscillations  are 
observed  in  the  free-stream  region  in  the  bottom  half  of  the  test 
section  that  are  nearly  in  phase  with  those  in  upper  free  stream. 

Measurements  were  made  over  multiple  planes  as  shown  in 
Figure  2.  The  streamwise  planes  yielded  the  (u,  v)  velocities 
in  the  (x,  y)  plane  as  well  as  the  span  wise  vorticity,  The 
cross-stream  planes  yielded  the  (v,  w)  velocities  in  the  (y,  z) 
plane  as  well  as  the  streamwise  vorticity,  co^.  The  coordinate 
system  for  these  measurements  is  also  shown  in  this  figure. 
Approximately  600-800  simultaneous  velocity  measurements 
are  made  over  each  plane.  The  results  described  in  this  paper 
represent  only  a  small  portion  of  this  data  set  related  to  the 
issues  addressed  herein.  The  complete  details  of  the 
development  of  the  velocity  and  vorticity  field  of  this  flow  can 
be  found  in  Cohn  (1999). 

VELOCITYA/ORTICITY  MEASUREMENT  METHOD 

The  velocity  field  was  measured  using  Molecular  Tagging 
Velocimetry  (MTV).  MTV  is  a  full-field  optical  diagiostic 
which  allows  for  the  non-intrusive  measurement  jS^f  the 
velocity  field  in  a  flowing  medium.  This  technique  can  be 
thought  of  as  the  molecular  equivalent  of  Particle  Image 
Velocimetry  (PIV).  Rather  than  tracking  particles  placed  in 
the  flowing  medium,  molecular  regions  are  tracked  by 
recording  the  luminescence  of  the  tracer  molecules.  A 
complete  description  of  the  implementation  of  this  technique, 
and  the  parameters  necessary  for  an  optimal  experiment,  can  be 
found  in  Gendrich  and  Koochesfahani  (1996),  Gendrich  et  al 
(1997),  and  Koochesfahani  etaL  (1997). 


Based  on  the  quality  of  the  images  and  conditions  in  the 
current  experiments  and  the  work  of  Gendrich  and 
Koochesfahani  (1996),  it  is  expected  that  the  instantaneous 
displacement  of  tagged  regions  is  measured  with  a  95% 
confidence  interval  of  less  than  0.1  pixel.  This  corresponds  to 
uncertainties  of  0.15  cm/s  and  0.2  cm/s  for  the  instantaneous 
streamwise  and  cross-stream  plane  velocity  measurements, 
respectively. 

The  0)^  and  fields  are  calculated  using  a  2"*^  order  central 
(4^^  order  accurate)  finite  difference  technique  from  their 
respective  velocity  fields.  The  use  of  standard  uncertainty 
analysis  indicates  that  the  uncertainty  in  the  vorticity  field  is 
given  by  6co  =  1.34  bxJh,  where  h  is  the  grid  spacing  used  in 
the  differentiation  and  6u  is  the  uncertainty  in  the  velocity 
measurement.  Applying  the  parameters  from  the  current 
experiment,  the  95%  confidence  level  for  the  instantaneous 
vorticity  in  the  current  experiments  corresponds  to  1.6  s'*  for 
(o^and  1.8  s'‘  for 

Previous  studies  of  this  forced  wake  suggest  that  many 
features  of  the  flow  are  quite  repeatable  and  phase-locked  to 
the  forcing  frequency.  This  allows  a  phase-averaging 
procedure  to  be  used  to  relate  data  from  different  fields  of 
view.  Results  presented  in  Cohn  (1999)  show  that  the  phase- 
averaged  data  is  a  very  good  representation  of  the 
instantaneous  result  in  this  flow.  Added  benefits  of  phase¬ 
averaging  are  the  reduction  of  the  measurement  uncertainty 
and  an  increase  the  “effective”  temporal  sampling  rate  of  the 
measurements. 

Using  the  bellow  command  signal  as  a  reference,  the  velocity 
field  data  are  phase  ordered  and  placed  into  64  non¬ 
overlapping  phase  bins.  For  each  field  of  view,  1000 
consecutive  velocity  realizations  are  acquired  at  a  sampling 
rate  of  33.3  ms.  There  are  approximately  16  instantaneous 
realizations  in  each  bin,  which  are  then  averaged.  With  1 6 
realizations  per  bin,  the  measurement  uncertainty  is  reduced  by 
a  factor  of  4  compared  to  the  instantaneous  values  mentioned 
earlier.  Using  64  bins  per  forcing  cycle,  the  effective  temporal 
spacing  between  consecutive  phase  bins  is  equivalent  to  2.6  ms 
in  real-time. 


SIDE  VIEW  RESULTS  AT  CENTER  SPAN 

Figure  3  illustrates  the  downstream  development  of  the 
phase-averaged  spanwise  vorticity  field,  <o)2>,  for  four  forcing 
amplitudes  as  weU  as  an  instantaneous  map  from  the  unforced 
case.  Note  that  each  forcing  case  combines  data  from  5  planes 
covering  different  streamwise  ranges  using  the  known  phase 
information.  The  vorticity  level  is  indicated  by  the  grayscale 
level,  and  also  by  contour  lines  placed  every  ±5  s'^ 

In  the  unforced  case,  the  vorticity  shed  from  the  splitter 
plate  turns  into  an  organized  vortex  array  at  about  x  =  6  cm. 
However,  the  horizontal  spacing  between  the  structures  is 
evolving  until  approximately  10  cm  downstream.  By  x  =  17 
cm,  the  magnitude  of  the  vorticity  within  the  spanwise  rollers 
has  decreased  to  a  very  small  magnitude. 

When  forcing  is  applied,  identifiable  vortical  structures  are 
apparent  at  the  earliest  upstream  location.  With  increasing 
forcing  amplitude,  the  peak  vorticity  level  of  the  vortical 
structure  increases  in  the  upstream  locations.  In  all  cases  the 
peak  vorticity  decreases  as  the  flow  moves  downstream  with  a 


rate  of  decrease  which  depends  upon  the  forcing  amplitude.  In 
fact,  at  the  farthest  downstream  locations,  the  peak  <g)2>  of  the 
high  amplitude  cases  are  less  than  that  of  the  low  amplitude 
cases.  This  is  seen  more  clearly  in  Figure  4  which  shows  the 
effect  of  downstream  location  and  forcing  amplitude  on  the 
peak  levels  of  <a)p>.  When  interpreting  the  absolute  values  of 
the  measured  vorticity,  it  should  be  noted  that  the  vorticity  is 
underestimated  due  to  the  spatial  averaging  inherent  in  the 
measurement.  This  mean  bias  error  has  been  estimated  in 
Cohn  (1999)  to  be  about  18%  for  the  conditions  of  these 
experiments. 

The  downstream  decrease  of  in  Figure  4b  is  not 
primarily  the  result  of  viscous  diffusion.  For  the  forcing 
amplitude  of  11.6%  the  estimate  of  the  effect  of  viscous 
diffusion  is  given  in  Figure  4b,  Also  given  is  the  estimate  of 
the  effect  of  vorticity  stretching  term  [(0w/0z)  caj  in  the 
vorticity  transport  equation.  A  value  of  dv//dz  =  -1  s‘^  was 
used,  which  is  consistent  with  our  measurements  of  this 
quantity  at  center-span.  Note  that  this  corresponds  to 
“negative”  stretching  and  a  reduction  of  the  vorticity  level.  It 
is  clear  in  Figure  4b  that  the  decay  of  peak  vorticity  by 
negative  stretching  dominates  over  viscous  diffusion  and  is 
consistent  with  the  measured  decay. 

The  measured  u(x,y),  v(x,y)  data  can  be  used  to  estimate  the 
value  dw/0z  based  on  the  continuity  equation  in  order  to 
provide  an  indication  of  the  extent  of  three-dimensionality  in 
the  flow.  The  uncertainty  in  the  instantaneous  value  of  dw/0z 
is  the  same  as  that  for  the  vorticity  described  earlier.  The 
downstream  variation  of  the  mean  6w/0z  at  center-span  is 
depicted  in  Figure  5  for  the  case  of  11.6%  forcing  amplitude. 
Comparing  this  result  with  the  earlier  mixing  data  suggests  that 
the  increase  of  three-dimensionality  indicated  by  the  increase 
in  the  magnitude  of  dw/dz  is  consistent  with  mixing 
enhancement.  Note  that  mixed  fluid  fraction,  begins  to 
increase  just  after  the  magnitude  of  the  mean  dv//dz  increases. 


END  VIEW  RESULTS 

Figure  6  compares  the  mean  streamwise  vorticity  field  for  the 
unforced  case  as  well  as  for  the  7.3%  and  11.6%  forcing 
amplitudes  at  four  different  downstream  locations.  Only  the 
right  half-plane  of  the  test  section  is  shown.  In  the  unforced 
case,  the  instantaneous,  and  thus  also  the  mean  values  are 
very  small.  For  the  7.3%  forcing  amplitude,  a  pair  of  counter¬ 
rotating  regions  of  vorticity  are  seen  near  the  sidewall  located 
at  z  =  -4  cm.  Note  that  these  regions  remain  close  to  the 
sidewall  as  the  flow  moves  downstream,  and  there  is  no 
evidence  of  the  presence  of  co,  at  center-span.  Therefore,  the 
modest  molecular  mixing  enhancement  at  center-span  at  this 
amplitude,  see  Figure  1,  is  not  due  to  small  scales  generated  by 
an  increase  in  three-dimensionality.  Rather,  it  is  due  to  the 
additional  “jelly-roll”  mixing  interface,  as  seen  in  MacKinnon 
and  Koochesfahani  (1997),  which  is  caused  by  the  increased 
spanwise  vorticity. 

The  11.6%  forcing  case  also  initially  exhibits  a  pair  of 
counter-rotating  vortices  near  the  sidewall  with  a  much  higher 
vorticity  level.  However,  farther  downstream  multiple  regions 
of  0)^  are  apparent  and  they  move  towards  the  center  of  the 
test-section.  Both  of  these  effects  are  responsible  for  the 


increased  three-dimensionality  and  the  observed  mixing 
enhancement  at  these  high  amplitudes.  We  draw  attention  to 
the  fact  that  the  regions  of  streamwise  vorticity  are  much  more 
complex  than  the  simple  vortex  pair  model  of  Roberts  (1985). 
The  presence  of  multiple  regions  of  has  been  argued  by 
Cohn  (1999)  to  be  the  result  of  the  reorientation  of  the 
spanwise  vortex  array  convecting  by  a  given  spatial  location. 


CONCLUSIONS 

The  evolution  of  the  vorticity  field  in  a  forced  low  Reynolds 
number  confined  wake  has  been  presented.  It  is  found  that  the 
peak  vorticity  of  the  spanwise  vortices  decays  with 
downstream  distance  primarily  due  to  negative  stretching.  The 
increase  in  the  magnitude  of  mean  0w/3z  appears  to  be  a  good 
indicator  of  mixing  enhancement  for  the  case  studied. 

For  a  medium  forcing  amplitude  with  modest  mixing 
enhancement,  there  is  little  evidence  of  streamwise  vorticity  in 
the  center-span.  The  mixing  increase  in  this  case  is  the  results 
of  additional  mixing  interface  generated  by  the  spanwise 
vorticity.  In  the  high  amplitude  forcing  case,  the  combination 
of  higher  streamwise  vorticity  levels  and  the  inward  movement 
of  the  regions  of  (x>j^  lead  to  the  previously  recorded  molecular 
mixing  enhancement  at  center-span. 
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Figure  1:  Streamwise  variation  of  mixed  fluid  fraction  at  mid-span  for  different  forcing  amplitudes  (Koochesfahani  and  Nelson, 
1997). 


(a)  (b) 

Figure  2:  MTV  measurement  planes,  (a)  Streamwise  (side  view)  measurement  planes,  (b)  Spanwise  (end  view)  measurement  planes. 
Note  that  the  origin  of  the  coordinate  system  is  at  the  center  of  the  test  section  at  the  tip  of  the  splitter  plate. 


Figure  4:  (a)  Effect  of  forcing  on  peak  vorticity  magnitude  at  several  downstream  locations,  (b)  Variation  of  vorticity  magnitude 
for  3  forcing  amplitudes  with  respect  to  streamwise  location.  The  solid  line  indicates  the  decay  expected  in  the  1 1.6%  forcing  case 
solely  due  to  negative  stretching  while  the  dashed  line  indicates  the  decay  solely  due  to  viscous  diffusion. 


